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We consider theories where the dark matter particle carries flavor quantum numbers, and has 
renormalizable contact interactions with the Standard Model fields. The phenomenology of this 
scenario depends sensitively on whether dark matter carries lepton flavor, quark flavor or its own 
internal flavor quantum numbers. We show that each of these possibilities is associated with a 
characteristic type of vertex, has different implications for direct detection experiments and gives 
rise to distinct collider signatures. We find that the region of parameter space where dark matter 
has the right abundance to be a thermal relic is in general within reach of current direct detection 
experiments. We focus on a class of models where dark matter carries tau flavor, and show that the 
collider signals of these models include events with four or more isolated leptons and missing energy. 
A full simulation of the signal and backgrounds, including detector effects, shows that in a significant 
part of parameter space these theories can be discovered above Standard Model backgrounds at the 
Large Hadron Collider. We also study the extent to which flavor and charge correlations among the 
final state leptons allows models of this type to be distinguished from theories where dark matter 
couples to leptons but does not carry flavor. 



I. INTRODUCTION 

It is now well established that about 80% of the matter 
in the universe is in fact dark matter, rather than visible 
matter [1]. However, the masses and interactions of 
the particles of which dark matter is composed are not 
known. One simple and well-motivated possibility is 
that dark matter is made up of particles with masses 
close to the weak scale that have weak scale annihi- 
lation cross section to Standard Model (SM) particles. 
Dark matter candidates with these properties neatly fit 
into the 'Weakly Interacting Massive Particle' (WIMP) 
paradigm, and therefore naturally tend to have the right 
relic abundance to explain observations. 

The matter fields (Q, U c , D c , L, E c ) of the SM each 
come in three copies, or flavors, that differ only in their 
masses. This reflects the fact that the Lagrangian of 
the SM possesses an approximate U(3) 5 flavor symmetry 
acting on the matter fields, which is explicitly broken by 
the Yukawa couplings that generate the quark and lepton 
masses. An interesting possibility is that the dark matter 
field, which we label by %, also carries flavor quantum 
numbers, with the physical dark matter particle being 
the lightest of three copies. Several specific dark matter 
candidates of this type have been studied extensively in 
the literature, including sneutrino dark matter [2-7] in 
the Minimal Supersymmetric Standard Model (MSSM), 
and Kaluza-Klein (KK) neutrino dark matter [8] in extra 
dimensional models. Recent papers on flavored dark 
matter include [9-11], (see also [12, 13]). 

In this paper we first consider the general properties 



of theories where dark matter carries flavor quantum 
numbers and has renormalizable contact interactions 
with the SM fields, and study their implications for direct 
detection and collider experiments. We then focus on a 
specific class of models where dark matter carries tau 
flavor, and perform a careful detector level study of its 
prospects for discovery at the Large Hadron Collider 
(LHC). 

To incorporate three flavors of the dark matter field, 
the flavor symmetry of the SM is extended from U(3) 5 to 
U(3) 5 x U(3) x , if x is a complex field such as a complex 
scalar, Dirac fermion or complex vector boson. If instead 
X is a real field, such as a real scalar, Majorana fermion or 
real vector boson the flavor symmetry is extended from 
U(3) 5 to U(3) 5 x 0(3) x . The new flavor symmetry U(3) x 
(or 0(3) x ) may be exact, or it may be explicitly broken 
as in the SM. 

Our focus in this paper will be on theories where dark 
matter has renormalizable contact interactions with the 
SM fields. Consider first the case where these contact 
interactions include couplings to the SM matter fields. 
These must be of the form shown in Fig. 1(a). If the dark 
matter flavor symmetry is to be exact, the field <p that 
mediates this interaction must transform under U(3) x 
(or 0(3) x ). If this vertex is to respect the SM flavor 
symmetry, cj) must also transform under the SM flavor 
group. In such a scenario, the different flavor states in the 
dark matter multiplet are degenerate, and the observed 
dark matter in the universe will in general consist of all 
three flavors. 

Alternatively, this contact interaction, in analogy with 
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FIG. 1: Vertices that link dark matter to a) the SM 

the SM Yukawa couplings, could represent an explicit 
breaking of the flavor symmetry. It is this scenario that 
we will be primarily concerned with in this paper. In 
this case the simplest possibility is that the mediator <fr 
is a singlet under both the SXI and the dark matter flavor 
groups. Then, if the SM matter field that x couples to is a 
lcpton, there is an association between the different dark 
matter flavors and lepton flavors. Accordingly, we refer to 
this scenario as 'lepton flavored dark matter'. Similarly, 
we label the corresponding case where x couples to a 
quark as 'quark flavored dark matter'. In this scenario 
the fact that the SM flavor symmetries are not exact 
naturally results in a splitting of the states in the dark 
matter multiplet, the physical dark matter particle being 
identified with the lightest. 

A different class of theories involves models of flavored 
dark matter where direct couplings between \ and the 
SM matter fields at the renormalizable level are absent. 
Instead, the contact interactions of \ with SXI fields are 
either with the W and Z gauge bosons, or with the 
Higgs, and can naturally preserve both the SM flavor 
symmetries and the dark matter flavor symmetry. The 
general form of such vertices is shown in Fig. 1(b). 
Closely related to this are theories with interactions of 
the exactly same form, but where dark matter instead 
couples to a new scalar <j> or vector boson Z' , which then 
acts as a mediator between the SXI fcrmions and the dark 
sector. In such a framework, dark matter is in general 
not associated with either quark or lepton flavor. We 
therefore refer to this scenario as 'internal flavored dark 
matter'. 

Since the characteristic vertices of lepton flavored, 
quark flavored and internal flavored dark matter are 
distinct, their implications for phenomenology are very 
different. In the next section we consider each of these 
classes of theories in turn, and study their collider signals, 
as well as their implications for direct detection and flavor 
physics. We then focus on a specific model of tau flavored 
dark matter and show that its collider signals include 
events with four or more isolated leptons and missing 
energy that can allow these theories to be discovered at 
the LHC above SM backgrounds. We also study the 
extent to which flavor and charge correlations among 
the final state leptons allows models of this type to be 
distinguished from more conventional theories where the 
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matter fields and b) the SM gauge and Higgs fields. 



dark matter particle couples to leptons but does not carry 
flavor, such as neutralino dark matter in the MSSM. 

II. FLAVORED DARK MATTER 

A. Lepton Flavored Dark Matter 

We first consider the case where dark matter carries 
lepton flavor. The lepton sector of the SM has a U(3)l 
x U(3)e flavor symmetry, where U(3)l acts on the 
SU(2) doublet leptons and U(3)e on the singlets. This 
symmetry is explicitly broken down to U(I) 3 by the 
Yukawa interactions that give the charged leptons their 
masses. (We neglect the tiny neutrino masses, which 
also break the symmetry). The characteristic vertex of 
lepton flavored dark matter involves contact interactions 
between x and the SM leptons of the form shown in the 
Fig. 1(a). The corresponding terms in the Lagrangian 
take the schematic form 

X A a L A Xa 4> +h.c, (I) 

if dark matter couples to the SU(2) doublet leptons L of 
the SM or alternatively, 

Xjjfmt +h.C, (2) 

if dark matter couples to the SU(2) singlet leptons E c . 
Here A is a U(3)l flavor index while i is a U(3)e flavor 
index and a is a U(3) x flavor index. There may also 
be additional interactions between the dark matter fields 
and the SM of the form shown in Fig. 1 (b) , in particular 
when x transforms under the SU(2) gauge interactions of 
the SM. 

The particle <j> that mediates dark matter interac- 
tions with the charged leptons is necessarily electrically 
charged. If x is a fermion then cf) must be a boson 
and vice versa. Any symmetry that keeps x stable will 
carry over to <f), and so <j> cannot decay entirely into SM 
states, except perhaps on cosmological timescales if the 
symmetry is not exact but very weakly broken. The same 
symmetry ensures that lepton flavor violating processes 
involving the dark matter field, such as fi — > ej, only arise 
at loop level through diagrams such as the one shown in 
Fig. 2. For concreteness, in what follows we take x to be 
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FIG. 2: Potential contribution to /i — > ej from lepton 
flavored dark matter. 



a Dirac fermion and to be a complex scalar, and restrict 
our focus to the case where \ couples to the SU(2) singlet 
lepton field E c , as in Eq. 2. The generalization to the 
other cases is straightforward, and is left for future work. 



1. Flavor Structure 

In general the matrix A will contain both diagonal and 
off-diagonal elements, thereby giving rise to lepton flavor 
violation. The experimental bounds on such processes 
are satisfied if all the elements in A < 10~ 3 for ~ 200 
GeV, even in the absence of any special flavor structure. 
In spite of these small couplings such a theory can still 
lead to interesting collider signals, since can be pair 
produced through SM gauge interactions, and will emit 
charged leptons as it decays down to the dark matter 
particle. Unfortunately, however, couplings of this size 
are by themselves too small to generate the correct 
abundance for \, if it is to be a thermal relic. This 
is not necessarily a problem if \ transforms under the 
SU(2) gauge interactions of the SM, or more generally if 
the theory has additional vertices of the form shown in 
Fig. 1(b), since these other couplings can play a role in 
determining the relic abundance. However, if x is a SM 
singlet and has no sizable couplings beyond those in Eq. 
2, the elements in A must be of order unity to generate 
the observed amount of dark matter, and aligned with 
the lepton Yukawa couplings to avoid flavor bounds. 

The matrix A can naturally be aligned with the SM 
Yukawa couplings if this interaction preserves a larger 
subgroup of the SM flavor group than just overall lepton 
number. For example, if we identify the three flavors 
of dark matter with the electron, muon and tau flavors 
in the SM, alignment is obtained if A, and the dark 
matter mass matrix, respects the U(l) 3 symmetry of 
lepton sector of the SAL In other words, the U(3) x x 
U(3) 2 symmetry is explicitly broken by A, and by the 
SM Yukawa couplings, down to the diagonal U(l) 3 . This 
larger symmetry forbids lepton flavor violating processes. 

A more restrictive possibility is that the only source 
of flavor violation in the theory is the SM Yukawa 
matrix, which then constrains the coupling matrix A to be 
consistent with minimal flavor violation (AIFV) [14]. In 
this scenario, the dark matter flavor symmetry U(3) x is 
identified with either U(3)_e or U(3)l of the SM, and the 



matrix A respects these symmetries up to effects arising 
from the SM Yukawa couplings. 

If we write the lepton Yukawa couplings of the SM as 

yA l L A E1R +h.c, (3) 

then the Yukawa matrix tja 1 can be thought of as 
a spurion transforming as (3,3) under the SU(3)lX 
SU(3)b subgroup of U(3)lX U(3)#. Consider first the 
case where U(3) x is identified with U(3)e- Then 

\ a l X a E c i <P +h.c.^A i V^ +h.c (4) 

If the theory respects MFV the matrix A is restricted to 
be of the form 

A/ = (al+/3yty).\ ( 5 ) 

Here a and /3 are constants, and we are keeping only the 
first non-trivial term in an expansion in powers of the 
SM Yukawa couplings. We write the dark matter mass 
term schematically as 

Kb Q x«x /3 , (6) 

In this case MFV restricts m x to have the form 

[m x }i J = (m l + Amy^y). 3 . (7) 

where too and Ato are constants. Since the SM Yukawa 
couplings arc small, the various dark matter flavors have 
small splittings and couple in a flavor diagonal way 
with approximately equal strength to leptons of the SM. 
Either the tau flavored or the electron flavored state will 
be the lightest, depending on the sign of Am. 

We now turn to the case where U(3) x is identified with 
U(3) L . Then 

AaVEftf +h.c. ->AaV^ +h.c. (8) 

MFV restricts the matrix A to be of the form 

Xa = a VA l , (9) 

where again we are working only to the leading non- 
trivial order in an expansion in the SM Yukawa couplings. 
The dark matter mass term now takes the form 

[m x ] A B = (m t + Am yrf) / . (10) 

We see that in this case the three dark matter flavors 
are again close in mass, but their couplings to the SM 
fields, though still flavor-diagonal, are now hierarchical. 
In particular, if the relic abundance is determined by 
A, we expect that only the tau flavor can constitute 
dark matter, since the couplings of the other flavors are 
relatively small. 

In the two Higgs doublet extension of the SM where one 
doublet gives mass to the up- type quarks, and the other 
to the down-type quarks and leptons, the coupling matrix 
A and the dark matter mass matrix are constrained by 
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FIG. 3: Dark matter scattering off a nucleus through photon exchange. 



MFV to be of the same form as in the SM. Therefore the 
formulas above continue to apply. However, since the 
lepton Yukawa couplings can be significantly larger in 
two Higgs doublet models, the tau flavored dark matter 
state may be somewhat split from the electron and muon 
flavored states, which remain nearly degenerate. 



2. Relic Abundance 

If X is a thermal WIMP, its relic abundance is set by its 
annihilation rate to SM fields. If \ is a SM singlet, and 
its only interactions are those of Eq. 2, then the primary 
annihilation mode is through f-channel <f> exchange to 
two leptons. In the relevant parameter space, the matrix 
A is constrained by flavor bounds to be very nearly 
flavor diagonal, so that lightest state in the dark matter 
multiplet is associated with a specific lepton flavor. We 
assume that the splittings between the different states in 
this multiplet are large enough so that only the lightest 
state is stable on cosmological timescales, and constitutes 
all of dark mater. 

The relevant terms in the Lagrangian, written 
schematically in 4-component Dirac notation, take the 
form 



£ => \ [x(i 
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(11) 



Here x represents the physical dark matter state and 
I the corresponding lepton. We have suppressed flavor 
indices since the matrix A is constrained to be nearly 
diagonal in the relevant region of parameter space. Since 
the dark matter particle is non-rclativistic at freeze-out, 
annihilation is dominated by the lowest partial wave. In 
this limit 
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(12) 



where we have assumed that m x » me, so that the 
masses of the final state leptons can be neglected. 

If x is also charged under the SM SU(2) gauge inter- 
actions then new annihilation channels open up. Dark 
matter can annihilate into two W's, two Z's, and also 
into SM fcrmions through through s-channel Z exchange. 
We leave a study of this for future work. 



3. Direct Detection 

The direct detection signals of this class of theories 
depend on whether the dark matter particle x transforms 
non-trivially under the SM SU(2) gauge symmetry, or 
remains a SM singlet. If x is a SM singlet, the leading 
contribution to dark matter scattering off a nucleus arises 
from the loop diagrams involving leptons shown in Fig. 
3. 

In the region of parameter space of interest to current 
direct detection experiments, bounds on lepton flavor 
violating processes constrain the coupling matrix A to 
be flavor diagonal. Therefore the dark matter candidate 
carries the flavor of the lepton it couples to. In this limit 
the relevant terms in the Lagrangian are again those 
shown in Eq. 11. As explained in Appendix A, this 
coupling gives rise to three distinct types of interactions 
between dark matter and the nucleus, specifically a 
charge-charge coupling, a dipole-charge coupling, and a 
dipole-dipole coupling. 

The differential cross section for the charge-charge 
cross section is given by the expression, 
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where tun is the mass of the nucleus, v is the velocity of 
the dark matter particle and E r is the recoil energy of the 
nucleus. Note that this is a spin-independent interaction, 
and hence is enhanced by Z 1 the total charge of the 
nucleus. The form factor F(E r ) appearing here is the 
charge form factor of the nucleus. It has been measured 
explicitly to be in good agreement with the Helm form 
factor [16]. The coefficient b p is defined as 
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Here rat is the mass of the lepton in the loop, which 
has the same flavor as the dark matter particle. The 
leading logarithmic part of this expression was calculated 
in [17]. In the case of electron flavored dark matter, the 
mass of the lepton mi in Eq. 14 must be replaced by the 
momentum transfer \k\ in the process, which we take to 
be 10 MeV as a reference value. 
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FIG. 4: Direct detection and relic abundance constraints on lepton flavor dark matter for a)x e and b)x r , when 
to^ = 150 GeV. The area above the solid blue curve is ruled out by the new XenonlOO [15] data. The green dashed 
curves signify the parameters for which we obtain correct relic abundance. 



The magnetic dipole moment of the dark matter can 
also couple to the electric charge of the nucleus. This 
interaction is also spin-independent. 



do 
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1 
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F 2 (E r ) (15) 



Finally, the dark matter can couple to the nuclear 
magnetic dipole moment via a dipole-dipole coupling. 
This interaction is spin-dependent, and therefore does 
not get an enhancement for large nuclei. It takes the 
form 
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where S nuc is the spin of the nucleus, /J. nuc is its magnetic 
dipole moment, and Fr){E r ) is the dipole moment form 
factor for the nucleus. There are currently no explicit 
measurements of the magnetic dipole form factor. A 
discussion of various form factors and an approximate 
calculation can be found in [18] (and references therein). 
The magnetic dipole moment of the dark matter particle 
\x x is related to the model parameters by 

Note that there is also a potential charge-dipole contri- 
bution to the cross section, where the dark matter vector 
bilinear couples to the magnetic dipole moment of the 
nucleus, but this interaction is suppressed by additional 
powers of momentum transfer. 

The dipole-charge interaction is sub-dominant to the 
charge-charge interaction. The dipole-dipole coupling, 
being spin-dependent, is also sub-dominant. Conse- 
quently, we use the charge-charge cross section for placing 
limits. Then, 
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Here oo is the cross section at zero-momentum transfer, 
and \x is the reduced mass of the dark matter-nucleus 
system. 

The ratio A/m^ corresponding to a thermal WIMP is 
plotted in Fig. 4 as a function of the dark matter mass, for 
the tau flavored and electron flavored cases. The current 
limits from the XenonlOO experiment [15] are also shown. 
It is clear from the figure that the expected improvement 
in sensitivity of the experiment by an order of magnitude 
will bring a large part of the parameter space of these 
models within reach. 

If x does transform under SU(2), we expect that the 
leading contribution to the cross section for dark matter 
scattering off a nucleus will arise from tree-level exchange 
of the SM Z, provided x carries non-zero hypercharge. If 
X arises from a representation which does not transform 
under hypercharge, then it does not couple directly to the 
Z, and so this effect does not arise. In this scenario, loop 
diagrams involving W bosons generate a contribution to 
the cross section [19] that must be compared against 
the contribution from the lepton loop above in order to 
determine the leading effect. 



4- Collider Signals 

What are the characteristic collider signals associated 
with this class of theories? For concreteness, we limit 
ourselves to the case where x does not transform under 
the SM SU(2) gauge interactions, and is a SM singlet. 
Then the mediator <f> also does not transform under the 
SU(2) gauge symmetry. 

We focus on the scenario where dark matter couples 
flavor diagonally, and where the electron and muon 
flavored states in the dark matter multiplet are highly 
degenerate, as would be expected from MFV. In such 
a framework, the charged leptons that result from the 
decay of a muon flavored state to an electron flavored 
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FIG. 5: Signal topologies at colliders in models of tau flavored dark matter. 



one (or vice versa) are extremely soft, and would be 
challenging to detect in an LHC environment. For the 
purposes of the following discussion, we will assume that 
these leptons are not detected. However, the splitting 
between a tau flavored state and an electron or muon 
flavored one is assumed to be large enough that the 
corresponding leptons can indeed be detected. 

The mediators <f) can be pair-produced in colliders 
through an off-shell photon or Z. Each <f> can then either 
decay directly to the dark matter particle, or decay to 
one of the heavier particles in the dark matter multiplet 
which then cascades down to the dark matter particle. 
Then, if the dark matter particle carries tau flavor, the 
decay of each <j> results in either a single tau, or in two 
charged electrons or muons and a tau. Each event is 
therefore associated with exactly two taus, up to four 
additional charged leptons, and missing energy. These 
event topologies are shown in Fig. 5. If, on the other 
hand, the dark matter particle carries electron flavor, the 
decay of each <fi will result in either a solitary electron 
or muon, or in two taus and an electron or muon. We 
therefore expect two electrons, two muons or an electron 
and a muon in each event, along with missing energy and 
up to four taus. 



B. Quark Flavored Dark Matter 

Let us now consider the case where dark matter 
carries quantum numbers under quark flavor. The quark 
sector of the SM has a U(3)q x U(3) c/ x U(3) d flavor 
symmetry, where U(3)q acts on the SU(2) doublet quarks 
and U(3)[/ and U(3)jj on the up and down-type singlet 
quarks. This symmetry is explicitly broken down to U(l) 
baryon number by the SM Yukawa couplings. 

The characteristic vertex of quark flavored dark matter 
has the form shown in Fig. 1(a). The corresponding 
terms in the Lagrangian take the schematic form 



quarks U c , we have 



\A a Q A X a <P +h.c. 



(19) 



X a l X a U^ +h.c. 



(20) 



This is easily generalized to the case where dark matter 
transforms under U(3)r>. 



K a X a D c a ^ +h.c. 



(21) 



if dark matter couples to the SU(2) doublet quarks Q. 
Alternatively, if it couples to the SU(2) singlet up-type 



Here the index A represents a U(3)q flavor index while 
i is a U(3)(7 flavor index and a is a U(3)_d flavor index. 
The mediator <j> is now charged under both color and 
electromagnetism. For concreteness, in what follows 
we again take x to be a Dirac fermion and to be a 
complex scalar, and restrict our focus to the cases where 
X couples to the SU(2) singlet quarks U c or D c as in 
Eq. 20 and Eq. 21. The generalization to other cases is 
straightforward, and is left for future work. 



1. Flavor Structure 

Contributions to flavor violating processes, such as 
K — K mixing, arise at loop level through diagrams such 
as the one in Fig. 6. The experimental bounds on flavor 
violation are satisfied if all the elements in A < 10 -2 , for 
nitf, <~ 500 GeV. However, as in the lepton case, couplings 
of this size are by themselves too small to generate the 
correct abundance for x, if it is to be a thermal relic. 
This is not necessarily a problem if \ has additional 
interactions with the SM, since these may set the relic 
abundance. However, if x is a SM singlet and has no 
other sizable couplings, the elements in A must be of order 
unity to generate the observed amount of dark matter. 
In this case the interaction matrix A must be aligned with 
the SM Yukawa couplings if the flavor constraints are to 
be satisfied. 

For the matrix A to be naturally aligned with the SM 
Yukawa couplings this interaction must preserve, at least 
approximately, a larger subgroup of the SM flavor group 
than just baryon number. This constraint is satisfied 
if the couplings A are consistent with MFV. In this 
framework, the only sources of flavor violation are the 
SM Yukawa couplings, and the matrix A respects the SM 
flavor symmetries up to effects that arise from them. 
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FIG. 6: Potential contribution to K — K mixing from 
quark flavored dark matter. 



Consider first the case when dark matter couples to 
the up-type SU(2) singlet quarks U c as in Eq. 20. MFV 
can be realized if the dark matter flavor symmetry U(3) x 
is identified with one of U(3)u, U(3)q or U(3)d of the 
SM, and the matrix A respects these symmetries up to 
effects arising from the SM Yukawa couplings. As in the 
lepton case, we will work to the leading non-trivial order 
in an expansion in powers of the SM Yukawa couplings. 

The quark Yukawa couplings in the SM can be written 

as 



y A a Q A D c a H + y A >Q A Utrf + h.c. . 



(22) 



The up-type Yukawa matrix y can be thought of as 
a spurion transforming as (3,3,1) under the SU(3)qX 
SU(3) c/ x SU(3) D subgroup of U(3) Q x U(3),yx U(3) D , 
while the down-type matrix y can be thought of as a 
spurion transforming as (3,1,3). 

If U(3) x is identified with U(3)c/, MFV restricts the 
matrix A to be of the form 



V = (al + 0yty)/ . 
while the mass term for \ becomes 



TO, 



V - 



(toqI + Am y'y) 



(23) 



(24) 



As explained earlier, we are working to the leading 
non-trivial order in an expansion in powers of the SM 
Yukawa couplings. Since the top Yukawa coupling is 
large, the higher order terms are not guaranteed to 
be small. However, this will not qualitatively affect 
our conclusions. The crucial point is that the dark 
matter states that couple to the first two generations 
of SM quarks are nearly degenerate in mass, and the 
mixing between them and the state with top flavor is 
small, protecting against flavor violating processes. The 
physical dark matter particle is expected to be either up 
flavored or top flavored, depending on the sign of Am. 
If U(3) x is identified with U(3)q we have instead 

\ A l = ay A \ (25) 

The dark matter mass term now takes the form 



[m x ] a B = (m t + Am yy^ + Am yy^ 



(26) 



While the first two flavors of \ are again quasi-degenerate 
in mass, their couplings to the SM are now hierarchical 
rather than universal. It is the smallness of the SM 
Yukawa couplings of the first two generations and their 
small mixing with the third generation that protects 
against flavor changing processes. If the coupling A is 
to generate the correct relic abundance the physical dark 
matter particle must belong to the third, rather than the 
first generation. 

Finally, if U(3) x is identified with U(3)d we have 



and 



K l = a (tfy) a z 



KJa = {m t + Amtfy) a 



(27) 



(28) 



Once again the first two flavors are very close in mass, 
and their couplings to the SM hierarchical. In order 
to give rise to the observed density of dark matter the 
lightest particle must be bottom flavored. 

We now turn our attention to the case where dark 
matter couples to the SU(2) singlet down-type quarks 
D c as in Eq. 21. MFV can be realized if the dark matter 
flavor symmetry U(3) x is identified with one of U(3)d, 
U(3)q or U(3);7 of the SM. The corresponding formulas 
for the form of the coupling matrix A and the dark matter 
mass may be obtained by simply interchanging y and y 
in the equations above. If U(3) x is identified with U(3)d, 
the matrix A is constrained to be of the form 



\ a b ={at+(3y^y) a b 



while the dark matter mass term is now 



\ m x\a = («*ol + Am tfy)\ 



(29) 



(30) 



We see that the states corresponding to the first two 
generations are quasi-degenerate in mass and couple 
universally to the SM, while the third generation can be 
somewhat split. This fact, together with the small mixing 
between the third flavor of dark matter and the first two 
allows flavor constraints to be satisfied. We expect that 
either the bottom or down flavor will constitute dark 
matter. 

If U(3) x is instead identified with U(3)q, we have 



Aa q 



a Va 



(31) 



while the mass term is of the form 

[m x ] A B = (m l + Amyyt + Am yy^ \ (32) 

While the first two flavors are still nearly degenerate, 
the different flavors now couple to the SM hierarchically 
rather than universally. We expect that dark matter will 
be composed of third generation particles. 

Finally, if U(3) x is identified with U(3)[/ these formulae 
become 



Q 



(yH) 



(33) 
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and 



Xu,d 



u, d 



[ m xV = (rn l + Am y^y). 



(34) 



Once again the first two flavors are very close in mass, 
and their couplings to the SM hierarchical. In order 
to give rise to the observed density of dark matter the 
lightest particle must be top flavored. 



Xu,d 



U, d 



2. Relic Abundance 

If the primary interaction of dark matter with the SM 
is through Eq. 20 or Eq. 21, then the relic abundance is 
set by t-channel annihilation to quarks. The calculation 
in this case mirrors that of the lepton flavored dark 
matter. In cases where the dominant annihilation mode 
is kinematically forbidden (e.g. for the top quark), 
three-body final states or loop-suppressed processes may 
dominate. 

In the region of parameter space which gives rise to the 
observed relic abundance, constraints on flavor changing 
neutral current processes require that the interaction 
matrix A be closely aligned with the quark Yukawa 
couplings. We therefore limit our analysis to the case 
where A is consistent with MFV. Then, in the mass basis 
the dark matter candidate is associated with the flavor of 
the quark it couples to most strongly, and does not mix 
significantly with the other flavors. The relevant terms in 
the Lagrangian, written in Dirac 4-component notation, 
take the schematic form 

Cd^ [x(1 + 75)9 + g(i - 7s)x f ] • (35) 

Here x represents the physical dark matter particle and q 
the corresponding quark. MFV ensures that the coupling 
matrix A is flavor diagonal in the quark mass basis, 
allowing us to suppress flavor indices. In the limit that 
the masses of the final state quarks can be neglected, we 
find for the annihilation rate 



(av) 



327r(m2 +TO 2)2 



(36) 



The additional factor of 3 relative to the lepton case arises 
because of the three colors of quarks. 

If dark matter transforms under the SM SU(2) gauge 
interactions, other annihilation modes open up, and 
may play the dominant role in determining the relic 
abundance. We leave this possibility for future work. 



3. Direct Detection 

The direct detection signals of this class of models 
depend on the flavor of quark that the dark matter 
particle couples to, and on whether or not \ transforms 
under the SM SU(2) gauge symmetry. Consider first the 



FIG. 7: Diagram contributing to direct detection for 
dark matter coupling to first generation quarks. 



case where x IS a SM singlet. In the region of parameter 
space relevant to current direct detection experiments, 
the matrix A is constrained by flavor bounds to be closely 
aligned with the SM Yukawa couplings. We therefore 
concentrate on the case where A is consistent with MFV. 
The relevant terms in the Lagrangian are then again 
those in Eq. 35. 

MFV suggests that the lightest state in the dark matter 
multiplet carries either the flavor of a first generation 
quark, or a third generation quark. The direct detection 
signals are very different in the two cases. If dark matter 
carries up or down flavor, it can scatter off quarks in the 
nucleus at tree level by exchanging the mediator <j> as 
shown in Fig. 7. 

Starting from the interaction in Eq. 35 we can integrate 
out the field <j>, leading to the effective operator 



4m- 



-X(1 + 7 5 )9 9(1-7 5 )X- 



After Fierz rearrangement, this operator becomes, 
A2 -X7"(1-7 5 )X 97^(1 +7 5 )<Z- 



8m- 



(37) 



(38) 



The dominant contribution to direct detection comes 
from the spin-independent vector-vector coupling. The 
dark matter-nucleus cross section (at zero momentum 
transfer) in this case is given by [20] 



°0 = 4^\A + Zf 



00 



647rm^ 

M 2 A 4 
QAnm 4 , 



[2 A - ZY 



(39) 
(40) 



for dark matter coupling to up-type and down-type 
quarks respectively. Here \x represents the reduced mass 
of the dark matter-nucleus system, while Z and A arc 
the atomic number and mass number of the nucleus. For 
a given value of A, this cross section is much larger than 
in the leptonic case. In fact, as shown in Fig. 8, the 
region of parameter space where \ can be a thermal relic 
is already excluded by direct detection experiments. 

We now move on to the case where the dark matter 
carries the flavor quantum numbers of third generation 
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FIG. 8: Direct detection and relic abundance constraints on quark flavor dark matter for a) Xu and b) xt, when 
to^ = 150 GeV. The area above the solid blue curve is ruled out by the new XenonfOO data [15]. The green dashed 
curves signify the parameters for which we obtain correct relic abundance. 
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FIG. 9: Diagrams contributing to direct detection for dark matter coupling to third generation quarks. The 
scattering can be off a) gluons and b) quarks via photon exchange. As discussed in the text, the photon exchange 

dominates for the example considered. 



quarks. The contribution arising from Fig. 7 is now 
suppressed by mixing angles, and is expected to be sub- 
dominant. There is a possible contribution to the cross 
section arising from the one loop diagrams shown in 
Fig. 9(a), where x scatters off gluons in the nucleus. 
However, it turns out that this is also not a significant 
effect. To understand why, we again integrate out the 
mediator <j> at tree level to obtain the effective operator 
shown in Eq. 38. This operator allows dark matter to 
scatter off a pair of gluons through triangle diagrams 
involving quarks. In general both the vector and axial 
vector terms in Eq. 38 contribute to the cross section. 
However, the contribution from the vector term vanishes 



identically as a consequence of the charge conjugation 
symmetries of QCD and QED (Furry 's theorem) [21, 22]. 
The axial vector interaction couples dark matter to 
gluonic operators that are parity odd rather than parity 
even [21, 23]. The parity symmetry of QCD can be 
used to show that the matrix elements of these operators 
in the nucleus are either spin-dependent or velocity 
suppressed in the non-relativistic limit, and therefore do 
not contribute significantly to dark matter scattering. 

Therefore, the dominant contribution to the cross 
section arises from one loop diagrams of the same form 
as in the lepton case, but now with the quarks running in 
the loop (Fig. 9 (b)). The cross sections will be identical 
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except for factors of color and charge. As before, we only 
use the charge-charge interaction to calculate the bounds. 



o 



fz 2 



3AVQ 
647r 2 mi 



log ( - 



(41) 



where Q 



3' 



■A for top and bottom quarks respectively, 



and m q is the mass of the quark in the loop. We see 
from Fig. 8 that the cross section corresponding to a 
thermal WIMP is within reach of current direct detection 
experiments. 

When dark matter transforms non-trivially under the 
SM SU(2) symmetry, scattering processes via Z-boson 
exchange can give large direct detection signals. Con- 
sequently, these scenarios are expected to be severely 
constrained. However, if \ arises from a representation 
which does not transform under hypercharge, then it does 
not couple directly to the Z, and so this effect does not 
arise. Then the effects of loop diagrams involving W 
bosons [19] must be compared against the contribution 
above in order to determine the leading effect. 



4- Collider Signals 



These theories are closely related to those where new 
particles, such scalar boson <j> or vector boson Z' , 
that have couplings of exactly the same form as in 1(b), 
mediate interactions between the SM matter fields and 
the dark matter sector. One important difference is that 
these can potentially give rise to SM flavor violating 
effects, if their couplings to the SM fields are off-diagonal. 

In this scenario, the dark matter states corresponding 
to different flavors may be exactly degenerate, if the 
internal flavor symmetry is exact, or split, if the sym- 
metry is broken. The collider phenomenology is highly 
sensitive to both the splitting between states, and to the 
particles produced when heavier states decay to lighter 
ones. The heavier particles in the dark matter multiplet 
can be pair produced through their couplings to the Z, 
the Higgs, (j) or Z' , and can then decay down to the 
lightest state. The additional particles produced in these 
decays, if visible, together with missing energy, constitute 
the collider signatures. This is a natural framework for a 
hidden valley [24] where particles such as <f> or Z' are the 
portal to the hidden sector. In this scenario, decays may 
be slow on collider time scales, giving rise to displaced 
vertices, since the couplings involved can be small in a 
technically natural way. 



The collider signals of this class of theories differ 
depending on whether the dark matter particle couples 
primarily to the third generation quarks, or to the 
quarks of the first two generations. We will restrict 
our discussion to the case where the couplings of \ are 
consistent with MFV. Our results can be extended to 
the more general case without difficulty. The mediators 
4> can be pair-produced at the LHC through QCD, and 
each will decay down to the dark matter particle either 
directly, or through a cascade. If the dark matter particle 
belongs to the third generation, at the partonic level each 
event is associated with two heavy flavor quarks and up 
to four light quarks. If, on the other hand, it belongs 
to the first generation, we expect between zero and four 
heavy flavor quarks in each event, along with two light 
quarks. 



C. Dark Matter with Internal Flavor 

Finally we consider the possibility that dark matter 
carries a new internal flavor quantum number that is 
distinct from either quark or lepton flavor, and does not 
couple directly to the SM matter fields at the renormal- 
izable level. In this framework, the only possible direct 
interactions of \ with the SM fields at the renormalizable 
level are to the weak gauge bosons or to the Higgs as 
shown in Fig. 1(b). These interactions do not generate 
large new sources of quark or lepton flavor violation. The 
direct detection signals are very similar to those of the 
corresponding theory where dark matter does not carry 
flavor. 



III. COLLIDER SIGNALS OF TAU FLAVORED 
DARK MATTER 

In this section we discuss in detail the collider signals of 
a specific model in which the dark matter particle carries 
quantum numbers under tau flavor. For concreteness we 
assume that the dark matter particle is a Dirac fermion 
which is a singlet under weak interactions, and therefore 
does not transform under the SU(2) gauge interactions 
of the SM. 

Dark matter couples directly to the SM through 
interactions of the form 

C = E i X J <t> + h.c] . (42) 

i— e,/^,r 

where <p is the mediator, and Xe,n,r are the dark matter 
and its copies. This interaction fixes the SM quantum 
numbers of <f>, which is charged under the photon and 
the Z, but does not couple to the W. 

For concreteness, we consider two benchmark spectra 
that are consistent with MFV, with \ transforming under 
U(3)e- Then Xe and x M are expected to be nearly 
degenerate since the corresponding SM Yukawa couplings 
are very small. We assume that Xt is lighter than Xe or 
Xju, and constitutes dark matter. 

We label the first benchmark spectrum tFDMI, 

m XjE = 110 GeV (43) 
m XiM = 110 GeV 
m X T = 90 GeV 
m = 160 GeV 
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FIG. 10: Pair-production cross section for the mediator 
<f> at the 14 TeV LHC run as a function of its mass. 

The second benchmark spectrum we study has a lighter 
mediator, and therefore leads to a larger production cross 
section (see Fig 10). We label this benchmark spectrum 
rFDM2: 

m x , e = 90 GeV (44) 
m XiM = 90 GeV 
m X T = 70 GeV 
m<p = 150 GeV 

In these simple models, only the mediator <fi carries 
SM gauge quantum numbers, so dark matter events 
at colliders must arise from 4> + 4>~ production. The 4> 
particles then decay, either directly or via cascade decays, 
into SM charged leptons and the dark matter particle. 
Therefore, the characteristic signature of this model is 
leptons+MET. 

As discussed earlier, MFV restricts the matrix A to be 
approximately proportional to identity. Consequently, 
we take the couplings of different dark matter flavors 
to SM to be equal, their common value set by the relic 
abundance requirement. Collider signals are insensitive 
to this value. 

Since 4> + 4>~ production proceeds through Drell-Yan 
and <j) is a scalar, the <fi pair comes out in a p- wave, leading 
to a small cross section, on the order of 10 fb at the 
14 TeV LHC run. Therefore, we do not expect the early 
LHC data to be able to probe this model. In order to 
obtain reasonable signal over background discrimination, 
tens of inverse fb of data will be required. 



A. Signal topologies 

Signal events come in three distinct topologies (see 
Fig. 5). Each (j> can decay directly into the dark 
matter particle and a r, corresponding to a short chain. 
Alternatively, it can decay to one of the heavier particles 
in the dark matter multiplet, which eventually cascades 
down to the dark matter particle, creating a long chain. 



Therefore each event can be categorized as comprising of 
short-short, short-long or long-long chains. 

Since r's are difficult to identify, we implicitly restrict 
ourselves to i = e, fj, final states in this section when 
we talk about leptons. The events with the long-long 
decay chain topology have four-lepton final states (not to 
mention a pair of r's), which have small SM backgrounds. 
When x is Majorana rather than Dirac, the short-long 
chain will also include a like-sign dilepton final state 
which is a very clean signal. The r's in the event 
could also decay leptonically, giving rise to additional 
leptons. However, these leptons are generally softer than 
the primary leptons. We focus on the long-long decay 
topology as the most promising channel. 

In order to simulate signal and background events we 
use the usrmod utility of MadGraph/MadEvent [25, 26], 
and we use BRIDGE [27] for the Xe,p, decays. Pythia [28] 
is used to simulate parton showers and hadronic physics, 
and PGS [29] with the default CMS parameter set is used 
to simulate detector effects. 



B. Backgrounds 

While four-lepton final states are rare in the SM, the 
signal cross section is also small so we carefully consider 
the three leading sources of backgrounds and devise cuts 
to reduce them as much as possible. 



1. (Z/7)W(Z/ 7 )W 

One of the dominant backgrounds is production of two 
opposite-sign, same-flavor lepton pairs from either on- 
shell or off-shell Z's and photons. Any missing energy in 
this background arises from mis-measurement of lepton 
momenta, which is small. For the following contributions 
to this background, we choose the following cuts: 

• Z — > £ + £~: This is the dominant component in this 
background, which we reduce by imposing a i?-veto 
(described in the next subsection) 

• Z — > t + t~ — > £ + £~: Even though the Z is on- 
shell in this process, the Z-veto is not effective 
due to the presence of neutrinos in the final state. 
This contribution is small due to the leptonic r 
branching ratios. The leptons arising from t decays 
are also softer, which we reduce by demanding the 
leptons to be energetic. 

• Z*/ 7 * -)• £+£--. While the off-shell production 
cross section is much smaller than on-shell produc- 
tion, this contribution is the main one that remains 
after the Z-veto and lepton energy cuts. We impose 
a missing energy cut to reduce this background 
component. 
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2. tt{Zh) 



(*) 



This background process, while it has a three-body 
final state, has a cross section comparable to the above 
process which is purely electroweak. When both tops 
decay leptonically and the (Z/7)M goes to leptons, the 
final state is 4£+jets+MET. The Z-veto reduces the on- 
shell Z production, and we also impose a dijct veto 
(described in the next subsection) in order to reduce this 
background, since signal events will typically not have 
any hard jets. 



3. WW(Z/i)W 

This process is qualitatively similar to the above 
process, but has a much smaller production cross section 
because it is purely electroweak. On the other hand, 
there are no additional hard jets in these events, so they 
escape the dijet veto. Consequently, events which escape 
the Z-veto can fake the four-lepton signal very well. 
Demanding the leptons to be energetic and imposing the 
missing energy cut helps reduce this background. 



4- Backgrounds with fakes 

There are also backgrounds arising from jets that are 
misidentified as leptons. We find that provided the 
fake rates arc of order 10~ 3 or less, the irreducible 
backgrounds described above are the dominant ones. 



C. Cuts 

We use the following cut flow in order to maximize 
signal over background: 

• Lcpton cuts - We demand events with at least four 
leptons each with pT > 7 GeV. At least two of 
these leptons are further required to have E > 50 
GeV. 

• Dijct veto - We discard events with two or more 
jets of pT > 30 GeV each. 

• Z veto - We veto events if the invariant mass of any 
Z-candidate (a pair of same-flavor and opposite- 
charge leptons) falls within 7 GeV of the Z mass. 
This is a tighter Z-veto than is usually used, but we 
find that the loss in signal efficiency is more than 
compensated for by the background reduction. 

• Missing energy - We require at least 20 GeV of 
missing energy in each signal event. Since most 
backgrounds with high MET have already been 
eliminated by the previous cuts in the cut flow, 
we find that a mild threshold such as 20 GeV is 
sufficient. 



Dataset 


Event rate after cuts at 100 fb 
Lepton cuts Jet cuts Z veto 


-i 

MET 


rFDMl 


46.73 


42.83 


38.41 


35.01 


tFDM2 


75.39 


69.30 


63.26 


57.04 


i+ri+r 


1617.94 


1582.42 


140.30 


13.32 


tti+r 


89.57 


19.45 


4.92 


4.70 


wwt+r 


14.70 


13.98 


2.51 


2.51 



TABLE I: Signal and SM background event rates for 
processes yielding 4-lepton final states after each set of 
cuts is progressively applied (note that £ = e,/i). All 
numbers are reported for the 14 TeV LHC run and 
include detector effects. 



D. Results 

The signal and background events of each type that 
survive these cuts are listed in Table I. These re- 
sults show that it is possible to discover the tFDM2 
benchmark above SM backgrounds at 5<r significance 
with about 20 ftr 1 of data at the 14 TeV LHC run. 
A higher luminosity (<~ 40 fb -1 ) would be needed in 
order to distinguish the rFDMl benchmark from the 
SM background. Note that while we have based this 
expectation on statistical uncertainties only, we have 
been conservative in many other aspects. In particular, a 
requirement that each event have at least one r candidate 
would virtually eliminate all remaining backgrounds 
while reducing the signal only moderately. Furthermore, 
one could do better than a pure counting experiment by 
taking into account the charge and flavor correlations 
present in the signal, which are different than the 
backgrounds in order to further increase sensitivity. We 
will indeed use this approach in the next section where we 
consider how the FDM model could be distinguished from 
more conventional DM models where the DM particle is 
a flavor singlet. 

While ATLAS [30] and CMS [31] have already per- 
formed searches in multilcpton final states, considering 
the low cross section of the FDM benchmark model, they 
are not yet expected to have exclusion level sensitivity to 
this scenario. 



IV. DISTINGUISHING rFDM 

Multi-lepton events with large missing energy are fairly 
common signals in theories with neutral stable particles 
and partners to the SM leptons, which include a variety 
of dark matter models. We would to like understand 
whether it is possible to distinguish at the LHC the 
model of tFDM that we studied in the previous section 
from models with similar signatures but where the dark 
matter does not carry flavor quantum numbers. Clearly, 
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FIG. 11: Flavor (red solid) and charge (blue dashed) correlations are shown for topologies in strawman models. 




FIG. 12: Flavor (red solid) and charge (blue dashed) correlations are shown for rFDM. Final state lepton charge 
ambiguities for Majorana dark matter models do not affect the charge correlation. 



this question is very difficult in general. Therefore, we 
focus on a more restricted question. We investigate 
whether it is possible to distinguish rFDM from a specific 
'strawman' model, where the dark matter does not carry 
flavor. 

The strawman model we choose is related to supersym- 
metric theories where the bino constitutes dark matter. 
The form of the lepton-slepton-bino vertex is very similar 
to the defining vertex of a theory of lepton FDM, except 
that in the supersymmetric case it is the slepton that 
carries flavor, not the bino. The strawman model we 
choose therefore consists of the bino, which we label by x, 
along with the three right-handed sleptons, Ef. The bino 
constitutes dark matter. To mimic the collider signals 
of rFDM, we add to the strawman model an additional 
'neutralino' which is heavier than the bino. x' is an 
admixture of a SM SU(2) doublet and singlet, so that 
it can be pair-produced through the Z, and is chosen to 
couple to leptons and sleptons in a flavor-blind way. This 
interaction takes the schematic form 

X'Ef X 'E ci + h.c. (45) 

The couplings of \' are somewhat different from those 
of a conventional neutralino in the MSSM, since any 
neutralino with significant couplings to the Z is expected 
to contain a significant Higgsino component, and the 
Higgsino does not couple universally to the different 
leptons. However, this simple strawman model captures 



the main features of theories where the dark matter does 
not carry flavor, while generating events which are very 
similar to those of rFDM. 

For simplicity, in what follows we assume that the three 
sleptons are degenerate in mass. In general, \' could 
either be lighter than or heavier than the sleptons, while 
the bino is the lightest of the new states. Both \' an d X 
are taken to be Majorana fermions as is the case in the 
MSSM. 

Signal events in the rFDM model involve four or 
more isolated leptons and missing energy. How does the 
strawman model generate similar events? The sleptons 
can be pair-produced in colliders. If they are heavier 
than x'i this leads to events of the form shown in 
Fig. 11(a), which involve six leptons, any or all of which 
could be taus. We label this possibility topology (a). 
Two x' particles can also be pair-produced, leading to 
events of the form shown in Fig. 11(b), which we label 
topology (b). These events involve four leptons, any or 
all of which could be taus. 

How can we distinguish between signal events in the 
two classes of models? One possibility is to note that 
we expect exactly two taus in each signal event in the 
rFDM model, whereas events in the strawman models 
will involve between zero and six. This could be a 
useful discriminant as the LHC experiments continue 
to improve their r identification capabilities. Presently, 
we do not make use of this discriminant. We also do 
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FIG. 13: The FDM spectrum and the strawman spectra compared 



not assume that the total event rate will be a reliable 
variable for discrimination. Even though in principle the 
event rate can vary widely across different models, for 
the following analysis we simply scale the event rate from 
the strawman model to match the number of events from 
tFDM and concentrate on ratios and asymmetries. 

In particular, we focus on charge and flavor correla- 
tions among the final state leptons in the event. In Fig. 11 
and Fig. 12, we exhibit the correlation of flavor and 
charge among the final state leptons in signal events for 
the tFDM model and the strawman model. The crucial 
observation is that, in the case of tFDM, for the chosen 
spectrum, the two upstream leptons arc also the hardest. 
This is likely to be the case for spectra motivated by 
MFV. These leptons arc charge anti-correlated since 
they arise from the decay of the charged mediators <f>. 
However, they have no flavor correlation, because the 
mediator does not carry flavor quantum numbers. 

Contrast this with the strawman model. Consider 
first events associated with topology (a). If the mass 
of <j> is much larger than that of \' , the two upstream 
leptons in the event are the hardest. These exhibit charge 
anti-correlation, but are flavor correlated, in contrast to 
tFDM. If, on the other hand, the mass of <fi is close to 
that of x' , t wo of the four downstream leptons will be 
the hardest. However, these exhibit no significant charge 
or flavor correlation, unlike rFDM. What about events 
associated with topology (b)? Here the two hardest 
leptons are again charge and flavor uncorrelated. We 
conclude from this that the charge and flavor correlations 
are different in the two theories, and may allow them to 
be distinguished. 

We generated signal events for the strawman model for 
three benchmark spectra, shown schematically in Fig. 13, 
and compared the resulting charge and flavor correlations 
to those of rFDM. In particular, the spectra we studied 
were the following. 



Spectrum 1 

We assume the masses of the sleptons to be the same 
as that of the mediator <p in rFDM. The masses of \' 
and x are & l so chosen equal to the Xe,^ and Xt mass 



respectively. Then, 



rn 



x' 



110 GeV 



(46) 



m x = 90 GeV 
iTiejiT — 160 GeV 



This spectrum can clearly give rise to both topologies 
in Fig. 11, but since the mass splitting between x & n d x' 
is small, events from topology (b) generally fail to pass 
the four-lepton cut requiring two leptons to have more 
than 50 GeV energy. Therefore, topology (a) dominates 
the phenomenology of this benchmark. 

In this topology, the two most upstream leptons are 
also the hardest, and are flavor-correlated. The tFDM 
leptons, as noted above, have no flavor correlation. 
Therefore, we expect that the flavor-correlation of the 
two hardest leptons is a good discriminant in this case. 
Spectrum 2 

If the mass of the sleptons is less than the mass of x', 
then only the topology (b) is allowed. The decay of x' is 
on-shell in this case. 

The representative spectrum we study is, 



160 GeV 
90 GeV 



me^T = HO GeV 



(47) 



In this case, the hardest leptons should exhibit neither 
charge nor flavor correlations, allowing us to distinguish 
it from rFDM. 
Spectrum 3 

Consider again the case when the mass of \' is less 
than the mass of the sleptons. As noted in the case of 
Spectrum 1, when the mass of x' is close to the mass of 
X, topology (a) dominates. On the other hand, when the 
mass of the x' IS very close to the mass of sleptons, the 
most upstream leptons become softer, and topology (b) 
dominates. The conclusions in this case are then identical 
to those of Spectrum 2. 

In the intermediate case, however, the result is a 
mixture of the two topologies. In order to investigate 
this we study a third spectrum, 



my = 140 GeV 
m x = 90 GeV 
me,ji,T = 160 GeV 



(48) 
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Datasct 



Frac. events with same 
Flavor Charge 



rFDMl 
tFDM2 
Spectrum 1(a) 
Spectrum 1(b) 
Spectrum 2 
Spectrum 3(a) 
Spectrum 3(b) 



0.52 
0.49 
0.87 
0.61 
0.55 
0.66 
0.60 



0.14 
0.14 
0.13 
0.39 
0.41 
0.33 
0.38 



TABLE II: Flavor and charge correlations for the two 
highest pT leptons in events passing cuts for different 
data samples. The strawman models are represented by 
the spectrum and their event topology. 



0.0 r 



-0.2 - 
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-0.6 - 



-0.8 - 



-1.0 L 



Spectrum 2 



Spectrum 1(b) 




Spectrum l(a 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 
(if 



In the next section we study the extent to which each of 
these spectra can be distinguished from rFDM. 

A. Comparison 

The correlations we obtain are listed in Table II. The 
results are in agreement with our expectations. Events 
with topology (a) in Spectrum 1 clearly exhibit flavor 
correlation between the two hardest leptons, as expected 
for the upstream leptons created from (flavor-carrying) 
sleptons. tFDM, on the other hand exhibits no flavor 
correlation in the hardest two leptons. 

In all the fake spectra with topology (b), the two 
hardest leptons show no preferential charge assignment 
beyond the ratio of 1 : 2 for same to opposite charge, as 
expected from random charge assignment. Consequently 
these cases have a weaker charge anti-correlation than 
the rFDM. 

Events from topology (a) in Spectrum 3 fall in the mid- 
dle, with somewhat significant charge anti-correlation, 
and a weak flavor correlation. While the correlation 
between charge and flavor is different from the rFDM 
case, higher statistics might be needed in this case to 
make a precise distinction. 

In Fig. 14, we plot the flavor and charge asymmetries 
of the two hardest leptons, along with statistical error 
contours for the rFDMl model at the 14 TeV LHC run, 
with 100 fb _1 of integrated luminosity. The asymmetries 
for flavor and charge are defined as, 



We see that with this amount of data, rFDM can be 
distinguished from the strawman models at 2a. 

Note that while we have relied only on the energies of 
the leptons to help us identify the topology of the event, 
this is just the simplest approach and can be extended 
with more sophisticated tools. For example, when 
applicable, the hemisphere algorithm [32] can distinguish 
particles arising from different decay chains, and has been 



FIG. 14: Flavor and charge asymmetry for different 
models and event topologies. We show la (green) and 
2cr (blue) contours for the rFDMl model at 100 fb _1 . 

The straight lines interpolate between points which 
correspond to different event topologies for each fake 

spectrum, in order to account for cases where both 
topologies contribute. 



used widely for this purpose [33-35]. There have also 
been other techniques proposed to help identify event 
topologies [36, 37]. And clearly, one can make use of 
further ratios involving the softer leptons to obtain better 
discrimination between these models. 



V. CONCLUSIONS 

In conclusion, we have studied the direct detection 
and collider prospects of theories where the dark matter 
particle carries flavor quantum numbers, and has renor- 
malizable contact interactions with the Standard Model 
fields. We have shown that the phenomenology of this 
scenario depends on whether dark matter carries lepton 
flavor, quark flavor or its own internal flavor quantum 
numbers. Each of these possibilities is associated with 
a characteristic type of vertex, leading to different pre- 
dictions for direct detection experiments and to distinct 
collider signatures. In particular, assuming a coupling 
consistent with relic abundance considerations, we have 
shown that many of these models could be probed in the 
near future by upcoming direct detection experiments 

We have studied in detail a class of models where 
dark matter carries tau flavor, where the collider signals 
include events with four or more isolated leptons and 
missing energy. We have performed a full simulation 
of the signal and SM backgrounds, including detector 
effects, and shown that in a significant part of the 
parameter space favored by MFV, these theories can 
be discovered above SM backgrounds at the 14 TeV 
LHC run. We have also shown that flavor and charge 



16 



correlations among the final state leptons may allow 
models of this type to be distinguished from simple 
theories where the dark matter particle couples to leptons 
but does not carry flavor. 



ACKNOWLEDGMENTS 

We thank Takashi Toma for pointing out an error 
in an earlier version of this manuscript. CK would 
like to thank Sunil Somalwar and Scott Thomas for 
helpful discussions about details of LHC searches and 
Rouven Essig, Jay Wacker and Eder Izaguirre for their 
generous offer to make their background event samples 
available to us. PA would like to thank Andrzej Buras 
for useful discussions and hospitality at the Technical 
University, Munich during the completion of a part of 
this work. We would also like to thank Michael Ratz 
for useful comments. PA and ZC are supported by NSF 
grant PHY-0801323 and PHY-0968854. SB acknowledges 
support from the CSIC grant JAE-DOC, as well as 
from MICINN, Spain, under contracts FPA2010-17747 
and Consolidcr-Ingenio CPAN CSD2007-00042. SB is 
also supported by the Comunidad de Madrid through 
Proyecto HEPHACOS ESP-1473. The work of CK is 
supported by DOE grant DE-FG02-96ER40959 and by 
the NSF Grant Number PHY-0969020. 



Appendix A: Direct detection of lepton flavored 
dark matter 



In this section we calculate the contribution to the 
cross section for dark matter scattering off a nucleus 
arising from the diagram shown in Fig. 3. 

Our approach will be to integrate out the mediator 
and the leptons / in the loop to obtain an effective vertex 
for the coupling of the dark matter particle \ to the 
photon. The resulting effective theory, which can be used 
to directly obtain the cross section, is valid provided the 
momentum transfer |fc| in the process is smaller than the 
mass of the lepton in the loop. The momentum transfer 
in direct detection experiments is typically of order 10 - 
50 MeV, which implies that this procedure is valid if the 
lepton in the loop is the muon or the tau, but not if it is 
the electron. However, the final result can generalized to 
obtain an expression that is approximately valid for this 
case as well. 

We first identify the operators that can potentially 
appear in the effective vertex. We begin by noting 
that the vertices in the diagrams in Fig. 15 do not by 
themselves violate CP symmetry. We therefore write 
down the leading effective operators that couple dark 
matter to the photon, and which are consistent with 
electromagnetic gauge invariance and CP. The lowest 
dimension operator consistent with these symmetries is 



unique. It is the dimension-5 dipole moment operator, 



X^X 



(Al) 



However, this operator does not actually appear in the 
effective theory. The underlying reason is that this 
operator breaks the chiral symmetry of the \ field. 
However, all the vertices and propagators in Fig. 15 
respect this symmetry, while the mass term of x, which 
breaks it, does not appear in the diagrams. This is 
most easily seen if we first integrate out just the heavy 
mediator <p an d consider the resultant effective four- 
fcrmion operator, 



X(1 + 75)«(1-75)X- 



(A2) 



A Fierz rearrangement shows that this is equivalent to 
the operator 



X7 A '(1-75)X^(1 + 75K, 



(A3) 



which establishes that the dark matter coupling is indeed 
chiral. Hence, we do not generate the dipole interaction 
above after integrating out the lepton. 

As a consequence, the leading contribution to the DM- 
nucleus scattering arises from dimension-6 operators. 
Again, gauge and CP symmetries, together with the 
chiral symmetry mentioned above, restrict us to the 
following two operators, 

01 = [x7"(l - r'Wx + h.c] F„„ (A4) 

2 = [ixVvl - 7 5 )^x + h.c] F°re^ p . (A5) 

The factors of i in the definitions of these operators 
have been chosen so that their coefficients in the effective 
theory are necessarily real. 

To calculate the coefficients of these operators in the 
effective theory we perform a matching calculation from 
the full theory to the effective theory, where the mediator 
<t> and the lepton I have been integrated out. 

In 4-component Dirac notation the relevant part of the 
Lagrangian in the full theory takes the form 

£ 3 \ [x(l + 75K <t> + *(1 ~ 7s)x f ] (A6) 

We can compute the one loop processes shown in Fig. 15 
to find the low energy effective Lagrangian. Since 
we are interested in direct detection processes with 
momentum transfer of at most 0(100 MeV), we only 
work to /m|) in momentum transfer. Further, we 
only keep the leading term in mtjm^ which is a good 
approximation in our case. In this limit, the amplitude 
in the full theory is given by, 



M 



X 2 e 



64n 2 m 2 



u(p 2 yfs{l - J 5 )u(pi)e*Jk) 



1 



■log 



-Api+P2) a kpe a ^ s 



(A7) 
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FIG. 15: Dark matter interaction with a photon in the full theory 




FIG. 16: Direct detection diagram for the dark matter 
in the effective theory 



where pi , p 2 and k are the momenta of the incoming dark 
matter, outgoing dark matter, and the photon, respec- 
tively Using integration by parts on the operators shown 
above, we see that the term with k 2 corresponds to 0i, 
and the second term corresponds to 02. Matching the 
coefficients, we can write down the effective Lagrangian, 



-eff 



64tt 2 



in: 



1 2 , 

2 + 3 l0g 



4 



(A8) 



We can now calculate the amplitudes for scattering of 
dark matter with nuclei arising from these different effec- 
tive operators and investigate their qualitative behavior. 
The scattering amplitude due to the first operator is 
given by, 



M 0l 



A 



2^2 



647T 2 77J 2 

q <t 



log 



x u( P2 )r(l -7 5 MPi) (N\Qq^q\N), (A9) 

where we sum the matrix elements of all quark bilinears 
in the nucleus, and Q is the charge of the quark in 
units of e. This is the typical interaction through the 
vector current. This gives rise to predominantly spin- 
independent cross sections which are enhanced for large 
nuclei. 

Consider the scattering amplitude due to the second 
operator, 



M 



o 2 



E 



iX e 



2 C 2 



327r 2 m 2 



u(p 2 h'"{l--f 5 )u(p 1 ) 



{P2+ / k f ka (N\Q^\N)e^ a0 . (A10) 



To disentangle different contributions, we use the Gordon 
identity on the dark matter spinors. Since we are using 
the equation of motion of the dark matter particle, we 
will now generate chiral symmetry- violating bilinears as 
well (ua a ^u in particular). Neglecting terms of higher 
order in momentum transfer and relative velocity, we get, 



M 



iX A e 



2 e 2 



(Pi +P2)u{P2+PlYk° 



647r 2 ra 2 i m x 4fc 2 
x u(p 2 ) [«cr" A1 7 5 ]w(pi) 

xJ2(N\Qqi f3 q\NK,ap ■ (Ail) 



We can rewrite a^^j 5 as ^crs p e"^ Sp and contract the 
Levi-Civita tensors. Using Gordon's identity again, the 
resulting expression can be brought to the following form, 



M c 



64tt 2 



in: 



■^2{N\Qqy a q\N) 



m x u(p 2 )<j af) u(p 1 )-j^ 



+ 2 U (P2)7 Q, u(pi) 



(A12) 



Combining both operators, the total scattering ampli- 
tude is 



M 



E 



ih' 



H x eu{p 2 )a ap u{ Pl )-^{N\Qq lm \N) 



+b p u{p 2 ) 1 f} u{p 1 ){N\Qq 1[j q\N)} 

where we have defined 

A 2 em v 



647r 2 m 2 
A 2 e 2 



b p = RA 2 2 (! + o lo § 
647r 2 mi 3 



) 



(A13) 

(A14) 
(A15) 



and neglected the velocity-suppressed contribution from 
M Ql - 

The first term in the amplitude corresponds to the 
magnetic dipole moment of x interacting with the nu- 
cleus, and the second term is the familiar charge-charge 
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interaction. The dipole couples to both the charge 
of the nucleus and its magnetic dipole moment. The 
momentum-transfer dependence of each of these terms 
is different. The dipole-charge interaction is enhanced 
at low-momentum transfers due to the presence of the 
k a /k 2 factor. However, the coupling to the dipole 
moment of the nucleon involves an additional power of 
the momentum transfer k. Therefore the dipole-dipole 
interaction has no such enhancement and exhibits the 
same recoil spectrum as the charge-charge interaction up 
to form factors. 

We show the three components of the scattering cross 
section: charge-charge (crzz), dipole-charge (<jdz) an d 
dipole-dipole {(tdd) [18, 38]. The differential scattering 
cross sections with respect to the recoil energy E r , are 
given as follows, 



dazz 2m 



dE r 



N 



4lTV* 



Z 2 bl F 2 (E r ) 



(A16) 



dapz _ e 2 Z 2 ^l 
dE r 4nE r 



1 - 



E r m x + 2m n 
v 2 2m N m x 



F 2 {E r ) (A17) 



dGDD 

dE r 



m N Mn«c M x 



S a 



3*SV> 



F 2 (E r ). (A18) 



Here tun is the mass of the nucleus, v is the velocity of 
the dark matter particle. S nuc is the spin of the nucleus, 
Unuc is the magnetic dipole moment of the nucleus, and 
Fd(Et) is the dipole moment form factor for the nucleus. 
The dipole-charge interaction is clearly enhanced at low 
momentum transfer. 



These results are only valid for the muon and the 
tau. However, in the case of the electron, the only 
significant difference is that it is the scale associated with 
the momentum transfer in the process \k\ that cuts off 
the logarithm in Eq. A15, and not the mass of the lepton 
m,£. In order to obtain approximate limits for the case of 
electron flavored dark matter it suffices to replace me in 
Eq. A15 by \k\. We choose = 10 MeV as a reference 
value. 
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